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Current reversals in a rocking ratchet: dynamical vs symmetry-breaking mechanisms 
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Directed transport in ratchets is determined by symmetry-breaking in a system out of equilibrium. 
A hallmark of rocking ratchets is current reversals: an increase in the rocking force changes the 
direction of the current. In this work for a bi-harmonically driven spatially symmetric rocking 
ratchet we show that a class of current reversal is precisely determined by symmetry-breaking, thus 
creating a link between dynamical and symmetry-breaking mechanisms. 
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Many processes in physics, chemistry and biology in- 
volve directed transport through periodic structures. For 
the equilibrium situation of Brownian motion, diffusion 
can be turned into directed diffusion by the applica- 
tion of a dc bias. In out-of-equilibrium systems, new 
mechanisms for directed transport may arise. Counter- 
intuitively, far from equilibrium it is possible to obtain 
directed transport through a macroscopically flat poten- 
tial in the absence of an applied dc bias. This is the 
so-called ratchet effect [iHll- 

The archetypal of a ratchet device is the rocking 
ratchet. In this setup, Brownian particles experience an 
asymmetric sawtooth potential and a sinusoidal rocking 
force. The rocking force drives the system out of equilib- 
rium, and directed transport is generated following the 
breaking of the symmetries of the system. An analogous 
effect can also be produced in a spatially symmetric po- 
tential and a biharmonic force, with the latter playing 
the double role of driving the system out of equilibrium 
and breaking the relevant time-symmetries |8l4l6|. 

A hallmark of rocking ratchets are current reversals. 
By progressively increasing the rocking force from zero, 
the generation of a current is observed, whose magni- 
tude is first an increasing function of the strength of the 
driving. However, at larger values of the rocking force 
the current reaches a maximum, then decreases to zero 
and changes sign. This feature can appear several times 
in a given system for different values of the force, thus 
producing multiple current reversals. Single and multi- 
ple current reversals have been observed in a variety of 
systems, both for an asymmetric potential and a sym- 
metric drive and for a symmetric potential and a time- 
asymmetric drive 18l-l20|. 



Current reversals are usually considered a dynamical 
effect, not related to the symmetry-breaking required to 
allow directed motion. In this work, for the specific sys- 



tem with a spatially symmetric potential and a time- 
asymmetric drive, we show that a class of current rever- 
sal is actually determined by dissipation-induced symme- 
try breaking. As a consequence, these reversals are not 
present in the Hamiltonian limit, nor in the overdamped 
limit. 

Our work consists into a theoretical analysis of the re- 
lationship between current reversals and dissipation in- 
duced symmetry-breaking. This is carried out compar- 
ing differing regimes: weakly damped, Hamiltonian, and 
overdamped. In the case of weak damping, where current 
reversals associated to dissipative effects are present, the 
theoretical analysis is also supported by experimental re- 
sults obtained with cold atom ratchets. 

In our theoretical analysis, the dynamics of particles 
in the considered spatially symmetric rocking ratchet is 
described by the Langevin equation 



U'{x) + F{t)+^{t), 



(1) 



where U{x) — Uq cos(2fca;)/2 is a periodic potential, a the 
friction coefHcient, ^(r) a Gaussian white noise: {£,{t)) = 
0, {^(t)^it')) = 2DS{t - t') und F{t) is an applied bi- 
harmonic drive of the form 241 



F{t) = Fq [Acos{ojt) + Bcos{2ujt + <j>)] 



(2) 



The generation of a current in such out-of-equilibrium 
set-up can be understood within the framework of the 
symmetry analysis ISMlSf. For the considered spa- 
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tially symmetric potential there are two time-symmetries 
which need to be broken to allow for the generation of a 
current: the shift symmetry, which corresponds to invari- 
ance under the transformation {x,p,t) — > (—x,—p,t + 
T/2), with T the period of the drive, and the time- 
reversal symmetry, which requires invariance under the 
transformation (x,p,t) — > (x, — p, —t). For a bi-harmonic 
drive of the form of Eq. [2l the shift symmetry is broken 
independently of the value of (j) {ioiA ^ 0, B ^ 0). The 
breaking of the time-reversal symmetry depends both on 
the value of the phase 4> and the dissipation level. In 



the Hamiltonian (dissipationless) case, the system is in- 
variant under time reversal for (f) = mr, with n integer. 
Thus, for these values of the relative phase, no current 
can be generated. However, for nonzero, weak dissipation 
the time-reversal symmetry is broken by dissipation, and 
directed motion can be produced also for cj) — mr. In 
the regime of weak dissipation, the dependence of the 
particles' velocity v on the phase (/) is well described, in 
leading order, by u = Wmax sin((/) — 0o) where ^o is deter- 
mined by dissipation, and vanishes in the Hamiltonian 
limit. It has been shown recently that this sine-like func- 
tional form is a consequence solely of the system symme- 
tries, being independent of the interaction details [21|. 
Of importance for the present study is the overdamped 
regime. In this limit, the so-called " supersymmetry" [15| 
(a;,p, t) ^ {x + A/2, —p, —t), with A the spatial period of 
the potential, is satisfied for </> = 7r/2 -I- m:, with n inte- 
ger. For these values of the driving phase no directed 
transport can occur. 

In order to establish a link between current reversals 
and symmetry breaking, we examine the dependence of 
the particles' current on two different quantities. First, 
we study the current as a function of the driving phase 
(/), which controls the time-symmetry of the Hamiltonian. 
This allows us to reveal the role of dissipation-induced 
symmetry breaking. Second, we consider the standard 
set-up for the observation of current-reversals: we fix the 
Hamiltonian by choosing a value of </> which corresponds 
to broken time-reversal symmetry {cf> = 7r/2, say) and 
study the current as a function of the driving strength. 
This will allow us to detect current reversals, and relate 
them to dissipation-induced symmetry breaking. In all 
the results presented in this work, the relative weight 
between the harmonic of the force is fixed {A = B ^ 
1 in all numerical simulations) and we vary the overall 
amplitude Fq. 

We examine first the weakly damped regime. The top 
panel of Figure [T] shows the average particles' velocity as 
a function of the driving phase for different values of the 
amplitude of the drive. The displacement of the current 
curves' extrema clearly indicates that, for a given dissi- 
pation level (i.e., a given value of the friction coefficient 
a), a variation of the drive amplitude leads to a varia- 
tion in the dissipation- induced phase lag (j)o- This is the 
central point of our analysis: because of dissipation, the 
curve of the average velocity vs the phase cj) acquires a 
nonzero phase lag (po, and the magnitude of the phase 
lag is a function of the strength of the driving. If we now 
examine the behavior of the current for a fixed phase (f> 
(say, = 7r/2) we see that the variation of (pQ due to the 
change in driving strength leads to a current reversal. 
Thus, a link between current reversals and dissipation- 
induced symmetry breaking is established. 

Our argument relies on the existence of a dissipation- 
induced symmetry breaking phase lag 0O: whose value 
depends on the dissipation level and driving strength. 
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FIG. 1: (Color online) Average atomic velocity as a func- 
tion of the relative phase (p between harmonics of the ac 
drive for several values of the driving amplitude Fq and 
oj = tJi, = k{2Uo/mY'^ . Top panel: simulation data in the 
weakly damped regime. The friction and the noise strength 
values are fixed to a = 0.15ao and D — 1.944Do, respectively, 
where ao = mkvo, -Do = ctQVo/k, and vo = (Uo/m)^'^ /lO. 
The values near (j> = 7r/2 (or (j) — 37r/2) show a current re- 
versal as the driving amplitude is increased. Bottom panel: 
simulation data in the Hamiltonian regime {a = D = 0). The 
initial conditions were chosen within the chaotic sea shown in 
Fig. [3l Lines are a guide to the eye. 



The consistency of the argument can be verified by con- 
sidering two extreme limits: the Hamiltonian case, and 
the overdamped regime. In both cases, the phase lag (po 
is "locked" to a given value (0 for the Hamiltonian case 
and 7r/2 in the overdamped regime) by the system sym- 
metry, and cannot be varied by modifying the driving 
strength. According to our argument, current reversals 
should disappear in both limits. 

We consider the dissipationless limit. The bottom 
panel of Figure [T] shows the average particles' velocity 
as a function of the phase cj) for different values of the 
driving strength. As in the weakly damped case, the 
current shows a sine-like dependence on the phase (p. 
However, unlike the weakly-damped case, the zeros of the 
current curves are now fixed by the time-reversal symme- 
try, and a change in the driving strength does not induce 
any phase shift. If we now consider the current depen- 
dence on the driving strength Fq for a given phase (say, 
(j) = 7r/2), we observe that no current reversal occurs. 
These results are consistent with the link we established 
between current reversals and dissipation-induced sym- 
metry breaking. 

It is interesting to study the dynamics of disappear- 
ance of current reversals while approaching the Hamil- 
tonian limit. This is done in Fig. [21 where for a fixed 
phase (p = tt/2 the current is studied as a function of 
the driving strength, for different levels of dissipation. It 
appears that by decreasing dissipation, the position of 
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FIG. 2: (Color online) Simulation data in the weakly damped 
regime: Current as a function of the driving amplitude Fo for 
several values of the friction a and a fixed driving phase (f} = 
7r/2. Rest of the parameters as in Fig. [l] As the dissipation 
is decreased, the current reversal position is shifted to lower 
values of the driving amplitude _fo- The diamonds show the 
results for the Hamiltonian system. Lines are a guide to the 
eye. 
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FIG. 3: (Color online) Simulation data in the Hamiltonian 
regime: The current as a function of the initial conditions 
w(0) and a::(0) is represented using a color density plot. The 
driving parameters are cj = ujv, Fq = O.lUok and (p — -k 12. 
A black circumference indicates the boundary of a regular 
circular region in which the current is zero throughout. 



the current reversal moves towards -Fq = 0) and in the 
Hamiltonian limit the current reversal disappears as its 
position coincides with Fg = 0. 

We notice that in the Hamiltonian limit the asymp- 
totic particles' velocity depends on the initial prepara- 
tion [17|. For the specific system of interest here, such a 
dependence is summarized in Fig. [3l There, we also ev- 
idence the chaotic sea. It is within that region that the 
initial condition for the calculations in the Hamiltonian 
limit presented here (bottom panel of Fig. [T] and one set 
in FiglS]) were chosen. 

We now consider the overdamped case. Numerical sim- 
ulations for this regime, in the presence of noise, are re- 
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FIG. 4: (Color online) Simulation data in the overdamped 
regime: Current as a function of the driving phase for sev- 
eral values of the driving amplitude Fo, with cj = 0.01cj„ and 
a friction a = lOOofQ. Top panel: D = 0.1944Do, being the 
lines a guide to the eye. Bottom panel: D — 0.1944 • 10~^-Do- 
The inset in the bottom panel shows the current as a function 
of Fq for a fixed driving phase </> = 1.01, displaying multiple 
current reversals. 



ported in the top panel of Fig. 2] The current is still 
described by a sine- like function v '^ sin(0 — (/)o) but the 
value of the phase lag is now locked by the symmetry 
to (^0 = 7r/2. A variation in the strength of the driving 
does not lead to a change in the phase lag. Studying 
the current as a function of the driving strength, for a 
fixed value (^ — 1:11 of the driving phase, reveals the ab- 
sence of current reversals. This confirms our statement 
about a link between current reversals and the symmetry 
breaking induced by the presence of a moderate amount 
of dissipation. 

So far we have demonstrated that the dependence 
of the dissipation-induced phase lag (/>o on the driving 
strength F^ results in current reversals, as observed by 
monitoring the dependence of the current at a fixed driv- 
ing phase on the strength of the drive. For the same 
argument, this type of current reversals are absent in the 
Hamiltonian and overdamped regime, where the phase 
lag 00 does not vary with the amplitude of the drive. It 
is important to specify the regime of applicability of such 
a reasoning, and underline that not all current reversals 
may be necessarily traced back to dissipation-induced 
symmetry breaking. In fact our argument relies on the 
assumption that the current can be well described by a 
smooth sine- like function v ~ sin((/)— (/iq). This holds pro- 
vided that the noise level is sufficiently large to smooth 
the curve. Thus, at low noise levels additional current 
reversals may appear, not related to dissipation-induced 
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FIG. 5: Experimental results for ID rocking ratchet for cold 
atoms, (a)-(f): average atomic velocity, rescaled by the recoil 
velocity Vr {vr = 5.88 mm/s for *^Rb), as a function of the rel- 
ative phase (j} between harmonics of the ac drive, for different 
values of the driving force amplitude, (g) Dissipation-induced 
phase lag 4>o, as obtained by fitting data as those in (a)-(f) 
with the function v/vr = Asm{cj> — (j)o), as a function of the 
driving amplitude. The amplitude A is reported in (h). The 
parameters of the optical lattice are: detuning from resonance 
A = — 15r and intensity per lattice beam II ~ 105 mW/cm^. 
The rocking force is of the form of Eq. (2] with u)/{2-n) = 100 
kHz, A= 1, B = 2 and Fq = —mfooj/k where m is the atomic 
mass, k the laser field wave-vector, and the value for fo (in 
MHz) for the different set of data is reported in the figures. 



symmetry breaking. An example of this is provided in 
the bottom panel of Fig. HI where current reversals are 
observed in the overdamped regime at a very low level of 
noise. 

The weakly damped regime, of central interest here, 
can be explored using dissipative cold atoms ratchets 
22|. We used a ID rocking ratchet set-up for ^^Rb atoms 
23j, which corresponds to a spatially symmetric poten- 
tial and a bi-harmonic rocking force of the form of Eq. [2] 
Proceeding along the lines of the theoretical analysis, we 
measured the atomic average velocity as a function of 
the driving phase 4>^ for a given dissipation level and for 
different strengths of the driving. Our results, shown in 



Fig. [5l confirm the prediction of the general theory we 
presented. An increase in the applied driving force am- 
plitude Fo leads to a large variation in the phase lag (po, 
as also summarized in Fig. ^g). Also the amplitude of 
the curve varies with the driving amplitude -Fo- How- 
ever, it never becomes zero (see Fig. [5jh)). It is thus the 
variation in the phase lag (jjo which produces a change in 
sign of the current observed for (p = ''i'/2, i.e. a current 
reversal. 

In conclusion, in this work we studied, both theoret- 
ically and experimentally, the relationship between cur- 
rent reversals and symmetry-breaking. For the specific 
system with a spatially symmetric potential and a time- 
asymmetric drive, we showed that a class of current rever- 
sal is actually determined by dissipation-induced symme- 
try breaking. As a consequence, these reversals are not 
present in the Hamiltonian limit, nor in the overdamped 
limit. 
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